Experimental Verification of the Quantized Conductance of Photonic Crystal Waveguides 
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We report experiments that demonstrate tiie quantization of the conductance of photonic crystal waveguides. 
To obtain a diffusive wave, we have added all the transmitted channels for all the incident angles. The conduc- 
tance steps have equal height and a width of one half the wavelength used. Detailed numerical results agree very 
well with the novel experimental results. 



Quantization of different physical quantities is one of the 
interesting phenomena in science. The quantization occurs 
because of the wave nature of particles. One quantity that 
gives strong quantization is the electrical conductance g, the 
inverse of the resistance. It is very difficult to calculate the 
transport properties of small devices analytically. Landauer 
was the first one to make the connection between the con- 
ductance and the transmission coefficient'. He showed that 
G - {e^ ITn:)T l{\-T), where T is the transmission coefficient. 
If one has a perfect metal T - \, then G ^ oo and therefore 
7? — > 0, as expected for a perfect metal. In 1981, Soukoulis 
and Economou^ proved, by using the Kubo-Greenwood for- 
mula, that for ID systems G - {e^/hn)T. There was a fierce 
controversy in 1980's regarding which formula was correct^, 
since the Economou-Soukoulis formula gives a finite value 
for the resistance for the perfect metal. The experiments^ re- 
solved this issue and indeed it was found that G = {e-/fin)T. 
Extensions to higher dimension^'^ were achieved for G and 
it was shown that for a multichannel wire the dimensionless 
conductance g 
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where f,„„ is the transmission coefficient between incident 
mode n and output mode m. The same set of modes A^^. are 
used for the incident and output modes. For classical waves 
the total T is equivalent to the dimensionless conductance^ in 
electronic systems T - g. For the ideal waveguide config- 
uration, g gives the number of propagation modes inside the 
waveguide. 

The key idea of Landauer was to relate the resistance of the 
sample with its transmission. So the idea of quantization of 
conductance can be also obeyed by all types of waves, elec- 
tromagnetic waves, acoustic and elastic waves. It is amazing 
that the only experiment'' that has been done with waves is 
the transmission of a slit of variable width for a given wave- 
length of A - \ .55 jim. Similar to its electronic counterpart, 
the optical conductance of a structure is described as the total 
light transmitted through the structure from a diffusive illumi- 
nation (an isotropic incoherent incident wave). So the con- 
ductance g for classical waves is dimensionless. In the exper- 
iment of Montie et al.'^ a two dimensional (2D) diffuser was 
used to achieve the diffusive illumination. The diffuser was 
essentially a 2D random array of scatterers through which the 



normally incident plane wave scatters diffusively and isotrop- 
ically. The diffused light passed through a metal slit and the 
transmitted light was collected. The result showed that the op- 
tical conductance increases in a staircase fashion. A new step 
occurs when the slit width W - nA/2 (n - 1,2,3, ...) i.e., a 
new mode is enabled in the slit. 

Photonic crystals (PCs) can be designed to have a bandgap 
that prohibits wave propagation in a certain frequency range'-. 
Line defects of PCs can confine light with frequencies within 
the bandgap inside the channel and act as waveguides. PC 
waveguides have been studied extensively and many appli- 
cations have been proposed'*'. However, the concept of op- 
tical conductance was not applied to PC waveguides until 
recently"'^ and no experimental work has been reported to 
our knowledge. In this letter we study the optical conduc- 
tance of photonic crystal (PC) waveguides, both numerically 
and experimentally. We will demonstrate experimentally that 
the optical conductance of a 2D PC waveguide in microwaves 
has the similar staircase effect as the metal slit. 

We design the working frequency of the waveguide at mi- 
crowave region, which makes the PC easy to be fabricated. 
One of the main experimental difficulties is the need of dif- 
fusive waves. It is easy to be achieved at optical wavelengths 
by a diffuser. For microwaves, a large enough random array of 
scatters are needed to act as a diffuser. Since the wavelength is 
now of the order of cm, the diffuser must be really big in size 
and the intensity of the transmitted waves will be too weak to 
be used efficiently. So we apply Eq. ([TJ in experiments for 
classical waves. We measure the transmitted power for each 
incident plane wave with different incident angle separately 
and sum them up to obtain the conductance of EM waves ' ' . 

The PC we study is a 2D square array of square alumina 
rods. The lattice constant is a = 11 mm and the square 
rods are of dimension d - 3.18 mm with relative permittiv- 
ity e = 9.8 and height h - \5 cm. It has a bandgap between 
9.43 GHz and 12.78 GHz for TM modes (Electric field par- 
allel to the rods). The waveguide is formed by two closely 
separated pieces of PC slabs with a size of 25fl x Aa each. The 
width of the waveguide (W) can be varied by changing the dis- 
tance between the two PC slabs. HP8510B network analyzer 
and a pair of horn antennas are used to measure the trans- 
mission. The antennas are mounted on motorized rotational 
stages which can move along a circle with a radius of 30 cm 
and centering at the middle of the entrance of the waveguide. 
So both the incident and outgoing angle can be controlled. 



-® 



Motor controlled 
linear stage 



Transmitter 



Motor controlled 
rotational stage 



Motor controlled 
yjptatlonal stage 



Motor controlled 
linear stage 



Receiver 



FIG. 1 : An illustration of the experimental setup. 
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FIG. 2: (Color online) The experimental and numerical results of 
optical conductance of a PC waveguide at 10 GHz. 



(See Fig. [T]for details). 

The experimental result of the optical conductance of the 
PC waveguide is shown together with numerical simulations 
in Fig. |2l We see that the optical conductance increases in a 
staircase manner with the increase of the waveguide width. 
The steps forming the staircases have essentially the same 
height and a new step appears at each integer multiple of 
A/2. Notice that the first conductance step doesn't start at 
W/A - 0.5. This is because that for the photonic crystal chan- 
nel the boundaries are not well defined and the field decays 
exponentially inside the photonic crystal. So the waveguide is 
wider than W. 

We use the commercial software COMSOL Multiphysics to 
do the numerical simulations. The simulation results of three 
frequencies are shown in Fig. [3] To understand the shape 
of optical conductance curves, the band structure of the PC 
with the channel are calculated using the supercell technique. 
When the channel width increases, some bands move down- 
wards into the bandgap of the PC. Fig.|4]shows the band struc- 
ture when W - 38.673 mm. There are four impurity bands 
inside the band gap. The experimental frequency / = 10 



FIG. 3: (Color online) The numerical results of optical conductance. 
The working frequencies are 10 GHz, 11 GHz and 12 GHz. The 
curves of 11 GHz and 12 GHz are shifted vertically. The PC has 4 
layers along the channel direction. Inset: g vs W/A when / = 10 GHz 
and the PC has 40 layers. 



GHz crosses the three impurity bands. It means the chan- 
nel can support three propagation modes. The propagation 
wave vectors k// of the modes can be read from the horizon- 
tal axis. Fig. |5] shows the relations between the propagation 
wave vectors of the propagation modes and the width of the 
channel when / = 10 GHz based on the band structure sim- 
ulations. The figure demonstrates clearly that the first propa- 
gation mode appears when W/A = 0.25 and one more prop- 
agation mode when the width increase A/2. Comparing with 
Fig. [3] it is clear that the optical conductance steps describes 
the number of propagation modes supported by the channel. 

The oscillation of the optical conductance curves can also 
be explained by the propagation wave vectors of the channel 
modes. A natural explanation about the conductance oscilla- 
tions is the Fabry-Perot interference, the interference between 
the multiple reflection of waves inside the waveguide. The 
phase difference between two succeeding reflection equals the 
propagation phase 2k//L (L is the length of the channel) plus 
the phase change of the reflections at the entrance and the exit 
of the channel. The reflection ratio changes with the channel 
width. To prove that the conductance oscillations are due to 
the Fabry-Perot interference, we calculated the optical con- 
ductance curves for a long channel (L - 4Qa, f = 10 GHz). 
The inset of Fig. [3] shows the simulation result. When the 
channel width changes from 0.41 A to 0.43 A, the oscillation 
has a full period. The reflection ratio is a continuous function 
of the channel width. So the reflection ratio will not change 
too much during the period since the width doesn't change too 
much. Then we neglect them and compare the propagation 
phases only. When W/A = 0.411,0.423,0.436, the conduc- 
tance curve has a dip, a peak and another dip. The propaga- 
tion phases are 2fc//L = 36.107r, 37. 097r, 38. 097r respectively. 
The phase differences between a dip and its neighboring peak 
are 0.99/: and 1 .OOtt. They agree with the Fabry-Perot expla- 
nation. 
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FIG. 4: (Color online) Band structure of PC with the channel. The 
channel width isW = 3S.673mm. The dark area is the band area of a 
perfect PC. The dash lines shows three working frequencies. The red 
solid curves shows the impurity bands in the band gap of a perfect 
PC. 
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FIG. 5: (Color online) Wave vectors of the propagation modes vs. the 
width of the channel. The solid lines are gotten from band structure 
of PC with the channel; The dash lines are extrapolation of the solid 
lines. 



From Fig. |2]we can see that the onset of the staircases in ex- 
periments and simulations are in good agreement. However, 
the Fabry-Perot oscillations are missing in the experiments. 
The major reason is that, while in simulations the model is 
perfectly 2D, we do not have an ideal 2D system in experi- 
ments. The source we get from a radiating horn antenna is a 
confined beam. When it encounters the PC structure, while 
the main part is propagating in xy plane, the beam also scat- 
ters out of the plane. Since the receiving antenna only collects 
the transmitted power in the xy plane where it sits, not all the 
transmitted power is obtained. At the onset of the staircase, 
more power is lost due to the multiple reflection. So the ex- 
perimental curve is lower than the one obtained by simulations 
and the Fabry-Perot resonances are missing. 



We also studied the optical conductance of a perfect electric 
conductor (PEC) channel, which serves as a simplified model 
of the photonic crystal waveguides. Analytical calculations 
starting from Maxwell's equations'^- are performed to under- 
stand the mechanism behind the stairs. 

Suppose y = is the interface of the air and the PEC chan- 
nel. The half-space y < is the air and y > is the PEC 
channel, x = is the middle of the PEC channel (see the inset 
of Fig. Hi. 

In TM mode, the incident wave is 



exp[i(k//y + kj_x)]. 



(2) 



Here k// - ko cos 0; k^ - ko sin 6. ko is the wave vector in air; 
6 is the incident angle. The reflected wave is the composition 
of plane waves. The total electric field in the air area is 
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E"" = exp[i{k//y + k^x)] + dkp(k) sxp[i(kx - kyy)]. (3) 



where ky - Jk^ - k^, and p(k) are unknown coeflicients. 

The channel is infinite long; so the field in the channel is 
the composition of all the channel modes moving towards +y 
direction. 



T^channel _ \ ^ a 



= X- 



„(x)exp[i^„>']. 



(4) 



Here 0m(x) - sin(m7r(xH-W/2)/W)is the channel mode profile 
andA„ = ^kl - {mnlWf. 

We have defined the fields before the channel (Eq. (O) and 
in the channel (Eq. (JUi) separately. £, and //j should be con- 
tinuous along the interface y - 0. By using these boundary 
conditions we can obtain the following set of equations for 
the unknowns A,,: 



Here 



^nynPn "^ ? .^m§mn ~ ^n- 



XW/2 
2k// exp[ik±x](l)„(x)dx. 
W/2 



XW/2 
<f)„(x)<p„(x)dx = W/2. 
w/2 



8mn=-z- I dk \ dx \ dx'lJk^^ - k^ 



(5) 



(6) 



(7) 



(p,n{x')^„{x) e,xp[ik{x - x')]\. 



(8) 



This method can be applied to TE mode too. Because of 
space consrains, we won't repeat it here. 
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FIG. 6: (Color online) Optical conductance of a infinite-long PEC 
channel. The first 10 channel modes are chosen in the analytic cal- 
culation; COMSOL Multiphysics is used to do the numerical simu- 
lations. 



The coefficients of the channel models A„ are solved by 
Equ. (|5]|. Then the optical conductance of the infinite long 
channel can be calculated. 

We also used COMSOL Multiphysics to simulate the infi- 
nite long channel numerically by inserting a perfect matched 
layer (PML)'"* in the channel. The incoming channel modes 
are absorbed by PML without reflection. Fig.|6]shows the ana- 
lytic calculation and numerical simulation results. The Fabry- 
Perot oscillations disappear since the channel is infinite long. 
The PEC channel has at least one TE propagaion mode no 
matter how narrow channel is. So the first step of TE mode 
begins at W/A = 0. It is an interesting question why the con- 
ductance curves, shown in Fig.|6], of TE and TM modes have 
different behavior 

Experimental measurements have also been performed and 
the results are shown together with the simulation results in 
Fig. |2]for TM mode. The experimental curve is smooth and 
lower than the simulation results, similar with the PC curves. 

We can simplify Equ. (|5]l furthermore. When koW ^> I and 
n <& koW, 
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FIG. 7: (Color online) The experimental and numerical results of 
optical conductance of a PEC channel for TM case. Inset: Diagram 
of the PEC channel. 



\jnPn if m^n; 
lO if m ?i «. 



(9) 



SoA„ X I„/2y„/3„. 

It means when the channel is wide, we can get the coeffi- 
cients of the first several modes directly. Then it is easy to 
prove that the contribution to the optical conductance by one 
mode converges to 1 when the channel width goes to infinity. 
This is the reason for the step profile. 

We have shown for the first time that the conductance of 
the photonic crystal waveguide is quantized in microwave fre- 
quencies. Quantization does not only occur in small length 
scales and in the quantum regime but also in long length scales 
and in the classical regime. We have also introduced new ways 
to obtain diffusive waves in the microwave region. 

This work was partially supported by Ames Laboratory 
(Contract No. DE-AC0207CH11385) and EU projects Meta- 
morphose and Phoremost. 
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